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Abstract 
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Q^ ' We study the power-law tails in the evolution of massless fields around a 

o ■ 
o 
l> 

' arguments for their existence at scri+, at the future horizon and at future 



fixed background geometry corresponding to a black hole. We give analytical 



timelike infinity. We confirm their existence with numerical integrations of 
the curved spacetime wave equation on the background of a Schwarzschild 
and a Reissner-Nordstrom black hole. These results are relevant to studies of 



^ ' mass infiation and the instability of Cauchy horizons. The analytic arguments 



also suggest the behavior of the full nonlinear dynamics, which we study 
numerically in a companion paper. 
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I. INTRODUCTION 



In the study of nonspherical gravitational collapse, the late stages of black hole formation 
and nonspherical stellar dynamics, certain simplifications have been traditional. Typically, 
linearized perturbation theory has been used on a fixed background, and the results have 
been taken to be representative of nonperturbative collapse. 

One of the most basic results about the evolution of fields on a curved background 
is that "sandwich waves" are not usually possible. At late times waves do not cut off sharply 
but die off in "tails." 

In the context of perturbations of spherical objects, like stars or black holes, arguments 
have been given leading to the conclusion that the tails have a specific power-law form. 
The intention of this paper is to analyze in detail in which regions of the spacetime this form 
for the tails holds and under what conditions they develop. 

In section II we start by giving an analytical outline of the development of tails in 
spherically symmetric fixed backgrounds. We present a somewhat more general and more 
pedagogical derivation of the results of the appendix of reference [||. Moreover, we give two 
new results: 

1) We show that for perturbations around a black hole, power-law tails develop not only 
at timelike infinity as was proven in [Q but also at scri^ and at the black-hole horizon. 
This result is of relevance since the development of tails in these regions is crucial for the 
physics of mass inflation |^ and the stability of Cauchy horizons. At least twice [§,0 in the 
literature it has been stated that the power-law tails are "nonradiative," suggesting that the 
tails made no appearance at scri^ or at the horizon. 

2) Generalizing the arguments for the Schwarzschild background, we show that power- 
law tails develop even when no horizon is present in the background. This would mean, 
among other things, that power-law tails should be present in perturbations of stars, or 
after the implosion and subsequent explosion of a massless field which does not result in 
black hole formation. 
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In sections III and IV we confirm the first of these results for the Schwarzschild back- 
ground by performing numerical integrations of the perturbation equations for different 
initial shapes and different multipole moments of the field. We also confirm that the results 
extend to the case of Reissner-Nordstrom black holes; this is the first clear evidence that 
power-law tails actually develop for the background of direct relevance to mass inflation 
and stability of Cauchy horizons. Section III contains a brief discussion of our numerical 
method; section IV contains our numerical results. 

In section V we make some final remarks, especially about possible implications, both 
for the behavior of test fields evolving on on a time-dependent stellar collapse or explosion 
background, and for the behavior of a (spherically symmetric) self-gravitating massless field 
in a collapse or explosion situation. We develop this subject in a subsequent paper. 

Tails have also been found important in the detailed calculation of gravitational wave- 
forms from the inspiral collapse of binary systems. The tail back reaction appear as a 
correction to the 3/2 Post Newtonian equations of motions 0. 

II. SCHWARZSCHILD BACKGROUND AND BEYOND 

In this section we examine the evolution of massless perturbation fields outside a star 
collapsing to form a black hole. We make the idealization that the mass loss in the transition 
from star to black hole can be neglected. (If there are massless perturbation fields outside the 
star at all, they will carry some energy away to null infinity, but we neglect this compared to 
the mass of the star.) If the star is spherically sjTumetric (again neglecting the perturbation 
fields, which need not be spherically symmetric), spacetime outside the star is Schwarzschild 
spacetime. We are therefore dealing with the (massless) wave equation on a region of 
Schwarzschild spacetime given by 

where we have used c = G = 1 units. We introduce a "tortoise" radial coordinate r^, = 
r + 2Mln(r — 2M), advanced time v = t + r^, and retarded time u = t — r^., in terms of 
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which the hne element becomes 

ds"" = -(1 - ^)dudv + r\de^ + sin^ ed<P^) . (2) 

We let the star begin its collapse at retarded time u = uq. The characteristic u = uq will be 
one boundary for our problem. 

After the surface has started collapsing, it rapidly approaches the speed of light. The 
wordline of the stellar surface is therefore asymptotic to an ingoing null geodesic v = vo- It 
is convenient to consider this v = vo null geodesic, rather than the stellar surface, as the left 
boundary of the problem, and the data on this line for u > uq to complete the specification of 
the problem. The variation of (p on the stellar surface is asymptotically infinitely redshifted. 
(See 1^ for details.) This means that (f),u will be small (more precisely, exponentially small) 
at late times and we make the specific assumption that, after some retarded time u = Ui, 
variations in on f = f o can be ignored. This black hole formation scenario is illustrated 
in Fig. 1. 

Since the background is spherically symmetric, each multipole of a perturbation 
field evolves separately. For a scalar field (f){t,r,6,ip), for example, we can write (p = 
J2i,m '^ini'^y ^)^™(^) f)/''^ and for each multipole moment we have 

- 4* . ^ V, . (l - ) (i(i±il + Rir)) . (3) 

where R{r) falls off as for large r, and where we have suppressed the indices l,m on 
The evolutions of scalar, electromagnetic, and gravitational perturbations are all governed 
by an equation with the form of (|^). Only the details of R{r) differ from one type of field 
to another. 

The general solution to (|^) can be written as a series depending on two arbitrary functions 
F and G: 

I oo 
p=0 p=0 

(4) 
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where the superindices on F and G indicate the number of times the function is differentiated; 
negative superindices are to be interpreted as integrals. The coefficients Af are dimensionless 
fixed numbers which have the same values as in flat spacetime. The coefficient functions 
Sf (r), however, vanish for M —>■ and have the form 



where the coefficients af are proportional to M. (See for additional details.) 

The most interesting characteristic data on u = Uq correspond to two different cases: a) 
an initially static perturbation field for u < uq, and b) a vanishing perturbation field for 
u < Uq. The first case corresponds to the collapse of an initially static star with an initially 
static perturbations. The second case corresponds to the collapse of an initially static star 
in which the collapse itself generates the perturbations. 

One may also consider a more general (nonstatic) initial perturbation on u = uq, but its 
generic effects can be modeled by the data on v = Vq, after refiection through the center. 
We shall see, on the other hand, that there is a qualitative difference in the tails if there is 
a static perturbation present initially. 

We break the problem of evolution of the fields into two steps. In the first step we find 
(f) on u = Ui. This step requires the data on -u = uq, but only the early (small u) data on 
the left boundary. In this step, then, the type of the left boundary (whether stellar center, 
surface of completely collapsing star, etc.) is unimportant. From the data on m = -Ui we find 
at late times, and at scri^, in a second step. To leading order in M, the evolution depends 
on the spacetime curvature in the first step, but not in the second step. The validity of this 
approximation is ultimately justified by numerical experiment. 

We begin with the first step, the scattering in the region uq < u < Ui. We have 
taken the variation in on f = fo be negligible after Ui. This means that, aside from 
backscattering, there is no outgoing radiation for u > ui. In (^) the first series represents the 
non-backscattered waves, so our assumptions about no late outgoing waves are equivalent 
to taking G{ui) = 0. For large r at u = Ui the dominant term in (^ is then 




(5) 
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^(m = Ml, r) = alr-^-^G^-^\ui) [1 + 0{M/r)] , (6) 

where G^~^\ui) = J^^ G{u)du. The characteristic data, then, is proportional to the inte- 
grated initial burst, as well as to M. In the case of an initial static field, with the form 

^static = /ir-'[l + 0(M/r)] , (7) 

the prediction is more definitive. In this case, at m = Mi, we have 

^ = M/x [(2/ + !)/(/ + 1)] r-('+i) [1 + 0(M/r)] . (8) 

(See again [H for details.) 

With this specification of characteristic data onu = ui we can go on to consider the sub- 
sequent evolution of fields. We will show that the late time behavior of the fields at constant 
r, and the late-time behavior at scri+, are independent of the details of the background at 
small r. The evolution is now the same in the spacetime geometry exterior to a star, a black 
hole, or to, e. g., an imploding-exploding shell of scalar field. 

We confine our attention to the region u > ui,r^ ^ M. The leading order effect on the 
propagation of \1/ is now that of the "centrifugal term" /(/ + l)/r^ in VJ. In this sense we 
now approximate spacetime as flat. The solution for \1/ is then that of @ with M = 0, 

I 

(To the approximation we are using here, for ^ M, we have also r = r^,.) By matching 
this form to the initial data on u = Ui we flnd that 

fiv) = (10) 

where 

Fo = {-iy2MG^-'\ui), P = 2, (11) 
if there is no initial static fleld,and 



(9) 
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Fo = (-2)'2M/x[/!/(2/)!], P = 1, (12) 

if there is an initial static field, with the form of (|^. (These expressions involve a sum over 
the Af, which are given in |^.) 

We next take t ^ r* and we expand g^''\u) = g^''\t) — g'^^^^\t)r^ + ■ ■ ■ , and similarly 
for /(f). By reordering terms we arrive at 

oo 

^=Y. K]^r:[f'-'''\t) + {-lYg^'^-\t)\. (13) 

n=—l 

The coefficients KJ^ are constructed from the A\ and vanish for / — n even, —l<n<l. We 
define h = f + {—i fg, and we note that for t ^ r* \E' has the form 

^ = Kr'Kty:' + Kr'^^h"{t)r:'+^ + ... (m) 



{t)r[ + Kf+^ \2 {t) - {t) 



We know that at late times f{v) has the power-law form given in (p!0[); we make the ansatz 
that g{u) also falls off as a power law, so that we may write h{t) = Ho/t^ , where Hq is a 
constant. 

We note that as t ^ oo the term h{t)r~^ will dominate if < P + 2/ + 1 and the term 
2j(2«+i)(^)^«+i ^jii dominate if iV > P + 2/ + 1. In either case only one of the constants Hq 
or Fq has an influence on the solution, and in fact the constant is simply an overall scaling of 
the solution. If either A^<P + 2/ + lorA^>P + 2/ + l, therefore, the form of the solution 
(aside from overall scaling) is fixed for ^ M, its continuation to smaller values of r* is 
fixed, and we cannot satisfy boundary conditions for small r (e.g., regularity at r = for a 
stellar model or, as r* — oo, at the horizon of a black hole). In order to satisfy a small-r 
boundary condition we must have = P + 2/ + 1 and hence 

f{t) ^ F,t-P, g{t) ~ (-l)('+i)Pot~^, (15) 



hit) ^ fit) + i-l)'git) ~ (16) 

For M -C r* <C t the form of the fields is particularly simple. In that case we have 
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^ = irf'+'2/(2'+i)(t)rl+i = -2Kf'+'Fo{P + 2/)!t-(^+2'+^Vl+i , (17) 

and \1/ therefore falls off as (initial static field) or (no initial static field) at 

timelike infinity i+. It follows from ([T5| ) and (^ that at scri+ (i. e., at f ^ oo) we have 

^(t; ^ oo, n) ~ A^g^^\u) ~ -A°Fo(/ + P - l)!?/"^-^ . (18) 

At null infinity scri+ the fields therefore fall off as (static initial field) or as (no 

static initial field). 

We keep in mind that the above analysis did not depend on the small-r details of the 
problem, and we go on to consider the specific case of a black hole. (It does not matter if 
it is eternal or formed in collapse, only that there is an internal "infinity", i. e. an event 
horizon.) As r^, — oo the curvature potential V/ in (|^) is negligibly small and the solution 
to d^) can be written as = a{u) + 7(f). The nature of the characteristic data at f = f 
requires that a{u) be a constant (aside from exponentially small variation) and with no loss 
of generality we choose it to be zero. For |r*| <^t we can then expand the solution, for large 
u and r=K <^ —M as 

^ = 7(y) = 7(t) + 7'(t)r, + ^Y{t)rl + ■■■ . (19) 

To join this solution, at r^, <^ — M to our previous solution in ([T7|) for ^ M , we make 
one further assumption. We assume that when t ^ |r^,| then, for the whole range of r* (from 

^ — M to > M), the solution has the form \E' ^ iinsta.t{r) /t^^^^^^ , where \E'finstat(''^) 
is a t- independent solution of (|^). This is clearly the case in the region t ^ ^ M. 
With this assumption we can match the solution in (|TP|) at r^, ^ —M and that in (|T^ for 

^ M, and conclude that 7(t) = Fot^^^^'^^, and therefore that the late time behavior at 
the horizon is 

^(m ^ 00, w) = 7(t;) = Tov'^-^^-^ . (20) 

The constant Fq is determined by the condition that there is a static solution \I'finstat such 
that 

8 



lim ^finstat(r) = To, (21) 

lim ^finstat(r) = -2Kf +^Fo(P + 2/)!r'+\ (22) 

The coefficient Fq, and therefore the behavior at the horizon, hke all other aspects of the 
late time behavior, is fixed by the initial backscattering. 

In the above analysis, we have seen that the backscattering of the initially outgoing waves, 
and the subsequent evolution in time, produces the power-law tails at scri^ and at future 
null infinity i+. The small-r nature of the background does not enter (except, of course, in 
the discussion of the tails at the horizon). This means that the same power-law tails should 
develop at scri^ and at i+ in models other than black hole models. We might consider as 
backgrounds: incompletely collapsing stars, static stars, imploding and exploding shells, and 
so forth. These different models would have different small-r boundary conditions. All that 
is important to tail formation is that the source of the perturbations is sharply cut off as 
happens, due to the infinite redshift, in the black hole collapse case. We might consider, as 
an example, a quadrupole deformation of a nonrotating neutron star. A dynamical process 
might change the quadrupole perturbation from one static value to another (nonzero) static 
value. In this case, clearly, there cannot be a power-law fall off of the perturbation in time; 
the perturbation does not fall off at late times. On the other hand there might be a phase 
change in which the stellar crust loses the ability to support shear stresses responsible for 
the quadrupole moment, and the star might quickly become spherical. In this case, our 
analysis predicts that at late times the exterior quadrupole perturbation will fall off as 

The only specific detail of our analysis that must be modified for non-hole models is the 
nature of the small-r solution. In particular, the horizon condition in ([19|), (pO|), and (pT]) 
must be replaced by the appropriate small-r condition, and the static solution ^finstat('") is 
no longer the solution well behaved at r — > 2M, but rather the solution satisfying the appro- 
priate small-r condition (e.g., regularity at r = 0). The result in (|2^ ) remains unchanged. 

The analysis of perturbations on a Reissner-Nordstrom (RN) background, of mass M 
and charge Q, is very similar to that given above for the Schwarzschild background. For 
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scalar perturbations or for all radiative degrees of freedom of the mixed electromagnetic- 
gravitational perturbations, the field equations can be put in the form of but with a 
potential of the form 



^ (r-r+)(r-r_) 



where r+ , r_ are the radii of the outer and inner horizons 



+ Rir) 



(23) 



r± = 2M(^l±^l-(g/M)2j , (24) 

and where R{r) falls off as r~^. The tortoise coordinate r* for the RN background is the 
solution of 

(25) 



dr (r — r^){r — r_) 
for r+ < r < oo. 

A review of the analysis of late time behavior for the Schwarzschild background confirms 
that almost all of the arguments depend on the general form of Viir) at large r and the 
exponentially sharp fall off of Viir) at the horizon. In these features there seems to be no 
difference between the Schwarzschild and the RN problems. The difference in the relationship 
of r and r* in the two cases must be carefully considered, however; it is this relationship [not 
the details of -R(r)] that determines the initial backscattering and therefore the behavior 
of the late time tails. It turns out that there is no difference (for the tails) between the 
Schwarzschild and the RN cases. When the expression in (4) is substituted in (3) (with the 
RN form of m, v and Vi) the result in (5) again emerges, and it is this result that determines 
the initial backscattering. The differences between the RN and Schwarzschild cases enter 
into the determination of B\{r) only to higher order in r~^. 



III. NUMERICAL METHOD 



It is straightforward to integrate equation (|^) on a uv null grid. The two-dimensional 
wave equation — 4\I' = V/(r)\E' is most simply discretized as 
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= + w^y - W5 - Axi V/ ( ^ 1 ^ hO^/i). (26) 

Here the points N, S, E and W form a null rectangle with relative positions indicated by the 
points of the compass, and h is an overall grid scale factor, so that 

Au = un — ue = uw — Us = 0{h), Av = vn — vw = ve — vs = 0{h). (27) 

Starting from null data on u = uq and v = vq, integration proceeds to the northeast 
(increasing v) on one u — const, line after another. 

The error estimate for follows directly if one assumes that the exact solution ^{u,v) 
has a Taylor expansion in the given null rectangle. On a grid of fixed total size in u and 
V there are 0{h~^) grid points, so that the total error on the far end of the grid from the 
null data should scale as h'^ when the grid size is scaled by h. We have done numerical tests 
which have confirmed this Oih?) convergence. 

The only difficulty in developing a code for this problem was the calculation of r from 

— {v — u)/2. We have done it by iteration of the defining equation, in the form r = 
r* — ln(r — 1) for large r, and in the form r\ = exp(r* — 1 — ri) for r close to 1, that is, close 
to the horizon, where ri = r — 1. We have used a grid of constant Au — Av. Because of the 
scale- invariance of the problem we have set 2M = 1. 

There is a small difference between the null data on u = uq and v = vo which we posed 
in the previous, analytical, section and our numerical work. For the numerical work, we 
posed constant data on v — Vq, and either generic (here Gaussian) or static data on u — Uq. 
The justification for this is the following. As we consider a linear wave equation, we can 
examine the evolution of data on v — Vq and of data on u — Uq separately. If we want to 
consider generic data on v — vq, i. e. data which are localized, we can refiect them back at 
the scattering potential to get generic (although different data) on some u = Uq. So if we 
are not very interested in the exact shape of our data, putting them on u — Uq is sufficient. 
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IV. NUMERICAL RESULTS 



We now report on the results of our numerical simulations of test fields on a Scliwarzscliild 
or RN background. Without loss of generality we have set 2M = 1. As we are dealing with a 
linear wave equation, the overall amplitude of our initial data can also be chosen arbitrarily 
and is physically irrelevant. As our initial data null surfaces we chose f = for u > and 
u = for t> > 0, which meet at the point t = 0, = 0. (The background is of course 
t-independent.) 

In a first series of numerical experiments, we chose Vi{r) appropriate for a massless 
minimally coupled scalar field on a Schwarzschild background for different multipole indices 
/. It is of the form @, with R{r) = 2M/r^. As null data we used a Gaussian of width 3 
centered at f = 10 on m = 0. \I' is chosen to be constant on the null boundary v = 0, with 
the constant determined by \I'(u = 0, w = 0). This is a simple approximation of the idea that 
the field is anchored in a star whose surface is rapidly redshifting, as explained in section II. 
We extended the grid from f = and u = to u = 400 and v = 400, with a typical value 
of Au = Av = 0.1. 

We examined the value of \E' as a function of t on three lines, namely r* = 10, m = 400 
and V = 400. We took these lines as finite approximations of the future timelike infinity 
of the future horizon H+ and future null infinity scri+. Log-log plots of these three "tails" 
for / = are shown in Fig. 2. The predicted power-law behaviors are apparent for the case 
of an initial static field. On and i+ the field falls off as t~^; on scri^ the fall-off is as 
The bend at the end of the scri+ line is not surprising. It represents null infinity only 
for V ^ u, which is no longer the case at that end of the line. The large wiggles at the left 
are a remnant of quasinormal ringing. 

Quasinormal ringing is shown more clearly in the linear plot, for / = 1 scalar pertur- 
bations, shown in Fig. 3. The agreement with the theoretically calculated frequency is 
good. From the plot we read off values of 0.56 and —0.19 for the real and imaginary parts of 
the frequency, for a scalar field with / = 1, for 2M = 1. The predicted values ^ are 0.58587 
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and —0.19532 for the least damped modes. 

In Fig. 4 we show the effect of the multipole index I on the power law of the tail at 
again for a scalar field and Gaussian null data. The agreement of the power laws with the 
prediction is striking. 

In a second series of numerical experiments we looked at the potential V/ appropriate for 
an electromagnetic field on a Schwarzschild background, which is @ with R{r) = 0. We 
took null data corresponding to the initial presence of a static /-pole moment. (These are 
given in power-series form in and we numerically evaluated the first few lowest orders 
in 1/r.) The tails at constant r for I = 1 and I = 2 are shown in figure 5. (There is no 
electromagnetic monopole field.) It was shown in how electromagnetic and spin-2 fields 
can be encoded in a scalar field. The potential V/ depends on the spin as well as on /, but the 
terms by which it differs for different spins are essentially Riemann tensor components and 
therefore of order 2M/r^. Terms of this order have been neglected in our analytic derivation, 
and the surprising prediction is therefore that the power law of the tails is independent of 
the spin of the testfield. 

In a third series of experiments we examined a scalar test field on Reissner-Nordstrom 
backgrounds. The potential for a scalar field is now of the form ([2B|), with R{r) = 2Mr~^ — 
2Q'^T~^ As explained in section III one expects the tails to be independent of the charge of 
the black hole. Figure 6 shows that this is in fact so, with the example of the tails at the 
(outer) horizon of an / = test field on Reissner-Nordstrom backgrounds with charge/mass 
ratios of Q^/M^ = 0.1 and 0.9 respectively. This figure constitutes also direct evidence 
for the existence of power-law tails on the (outer) horizon for a generic perturbations, thus 
underpinning a central condition for the mass inflation scenario. 

V. CONCLUSIONS 

When the dust of an approximately spherical collapse has settled, and spacetime inside 
the future lightcone of the collapse has approached Schwarzschild, Minkowski, or a stellar 
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interior spacetime, any massless fields that were present in the collapse still show "tails" 
that linger, decaying only as a power of time. In particular the l-th multipole moment of a 
massless test field decays like fixed radius at large times, with P = 1 if there 

is an Z-pole moment present before the collapse and P = 2 otherwise. But power-law tails 
are also present on scri^, where they decay like u~^''~^^\ and if a black hole has formed, on 
the horizon, where they decay like The amplitude of the tails can be calculated 

as a function of the initial multipole moment, and in its absence, as an integral over the 
radiation going out to infinity during the collapse. 

Neither this amplitude, nor the exponent of the power law, depend upon the spin of 
the massless field in question, nor, if the black hole is charged, on its charge. We have 
shown the origin of these features in an analysis which is based on [Q. The spin- and 
charge- dependent parts of the effective radial potential of a massless test field on a Reissner- 
Nordstrom background are of an order that can be neglected in this analysis. 

We have checked numerically that the tails are indeed present on null infinity and the 
horizon, and are independent of the spin of the field and the black hole charge. The fact 
that power-law tails are present on the outer horizon of a Reissner-Nordstrom black hole 
after a generic collapse situation is of crucial importance to the mass inflation scenario. It 
has, to our knowledge, never been demonstrated explicitly. 

Finally, one decisive step of our analysis was a regularity condition, either on the horizon 
when a black hole formed in the collapse, or else at the center of spherical symmetry. This 
argument generalizes to any kind of boundary condition posed at small radius, and strongly 
suggests that perturbations of massless fields outside any spherical background object should 
also have power-law tails with the powers given above. In a subsequent paper we report 
results for a closely related nonlinear problem: the implosion of a shell of scalar field. 
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FIGURES 



Fig. 1: Spacetime of a collapsing spherically symmetric star. Outside the star the metric 
is Schwarzschild. Coordinates t, r^, u and v are indicated. 

Fig. 2: Log-log plots of a spherically symmetric scalar test field on Schwarzschild space- 
time. The initial data were Gaussian. 0(r = 10, t) represents future timelike infinity. 
(I){u — 400, v) represents the future horizon. 0(v = 400, u) represents future null infinity. 
The power-law exponents are —3.08 on i+ and —3.02 the horizon, and —2.11 on scri+. 

Fig. 3: 0(r — 10, t) plotted for a scalar testfield with angular dependence Z = 2 on 
Schwarzschild spacetime, showing quasinormal ringing. The initial data are Gaussian. 

Fig. 4: Log-log plots of 0(r = 10, t) for scalar testfields with angular dependence / = 
0, 1, 2, showing power-law exponents —3.03. —4.99, and —6.94. Gaussian initial data. 

Fig. 5: Electromagnetic test field, with angular dependence Z = 1, 2, 3 on Schwarzschild 
background. Plotted is the field at it = 400 versus v, representing the horizon. Static-static 
initial data. The power-law exponents are —4.95 and —5.93 and —8.63. 

Fig. 6: Scalar test field on RN background, shown on "horizon". Gaussian initial data 
with I — 0. Shown are charge/mass ratios Q'^/M'^ — 0.1 and 0.9. The power laws are —3.18 
and -3.17. 
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